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ABSTRACT 


We  developed  a  Micro  SAINT  computational  program  that  executes  a  series  of  predictive 
equations,  developed  by  the  U.S.  Army  Research  Institute  of  Environmental  Medicine,  for  deep-body 
temperature,  heart  rate,  and  sweat  loss  responses  of  clothed  soldiers  performing  physical  work  at 
various  environmental  extremes.  The  user  can  employ  the  program  by  inserting  input  parameters  in 
the  scenario  section  of  the  MicroSAINT  program. 

The  developed  program  could  be  used  to  help  avoid  casualties  associated  with  environmental  heat 
extremes,  to  predict  appropriate  work-rest  cycles,  or  to  predict  water  requirements  due  to  sweat  loss. 
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1.  INTRODUCTION 


For  the  last  two  decades,  the  U.S.  Army  Research  Institute  of  Environmental  Medicine 
(USARIEM)  has  been  establishing  a  data  base  and  developing  a  series  of  predictive  equations  for  deep- 
body  temperature,  heart  rate,  and  sweat-loss  responses  of  clothed  soldiers  performing  physical  work  at 
various  environmental  extremes.  Predictive  equations  for  rectal  temperature,  heart  rate,  and  sweat  loss 
as  functions  of  physical  work  intensity,  environmental  conditions,  and  particular  clothing  ensembles 
have  been  published  in  the  open  literature  (1-11). 

In  this  report,  we  describe  a  Micro  SAINT  (12)  computer  model  based  on  the  USARIEM 
equations  for  deep-body  temperature,  heart  rate,  and  sweat  loss.  The  Micro  SAINT  program  is  written 
with  a  user  friendly  approach.  The  user  may  supply  the  input  parameters  in  the  scenario  section  of  the 
program  and  obtain  the  results  as  output  stored  in  Micro  SAINT  snapshot  files.  As  such,  the  user  may 
view  the  main  program  as  a  “black  box”  because  input  and  output  information  are  stored  separately  in 
the  scenario  and  snapshot  sections  of  the  Micro  SAINT  program. 

The  following  sections  constitute  a  detailed  summary  of  the  USARIEM  heat-stress  model,  the 
Micro  SAINT  model,  and  an  illustration  of  bow  to  use  the  Micro  SAINT  program.  The  Micro  SAINT 
program  is  listed  in  the  apnendix  and  can  be  obtained  from  the  authors  in  return  for  a  5  1/4  inch, 
double-sided,  double-density  disc  and  a  self- addressed  disc-mailing  envelope. 

2.  MATHEMATICAL  EQUATIONS 

EQUILIBRIUM  RECTAL  TEMPERATURE  AND  HEART  RATE  PREDICTIONS 
The  general  formula  for  predicting  the  final  equilibrium  rectal  temperature  (Trey) 
suggested  by  Givoni  and  Goldman  (1)  ..  •  '  Pandolf  et  al;  (6)  is: 

Tre/(°C)  =  36.75  +  0.004(M  -  Wcx)  +  0.0011  Hr+C  +  0.8  Exp[0.0047(Ereg  -  EJ]  (1) 

Equation  1  is  comprised  of  three  components:  metabolic,  dry  heat  exchange,  and 
evaporative  heat  exchange.  The  metabolic  component  is 

[37.75  +  0.004(M-Wei)] 

in  which 

M  =  1.5  W  +  2.0  (W+L)(L/W)2  +  i/(W+L)[1.5  V2W  +  0.350^]  (2) 

as  published  by  Pandoff  et  al.  (5),  and 

Wex  =  0.098G(W+L)VU,  (3) 

as  published  by  Givoni  and  Goldman  (1),  where 
M  =  metabolic  rate  (watts) 


Wex  =  external  work  (watts) 

W  =  nude  body  weight  (kg) 

L  =  clothing  and  equipment  weight  (kg) 

T]  =  terrain  factor  (defined  in  Table  1) 

=  walking  velocity  (m/s)  and 
G  =  grade;  (%) 

The  dry  heat  exchange  component  is 

[0.001 

where 

Hr+C  =  6.45  Ad(Tdb-Ttk)/lt 
as  given  by  Givoni  and  Goldman  (1) 

Aj  body  surface  area  (m2) 

Tjb  =  dry  bulb  temperature  (°C)  . 

Tifc  =  average  skin  temperature  (”C) 

It  =  total  insulation  including  air  layer  (Ia)  and  intrinsic  clothing  (IcJ) 

The  evaporative  heat  exhange  component  is 

0.8  exp  [0.0047  (Erc,  -  Emax)} 
as  given  by  Givoni  and  Goldman  [1]  where 

2rc,  =  (M-  WeI)  +  Hr+C 

and 

Emai  =  14.21  im/lt  kdcjj  (P.fc-daPa) 

where 

im  =  permeability  index  (N.D.) 

A deff  =  effective  surface  area  for  evaporation  (m2) 

=  water  vapor  pressure  at  the  skin  (mm  Hg) 

<j>a  =  relative  humidity  (%)  and 

Pa  =  saturated  water  vapor  pressure  of  air  at  Tdb  (mm  Hg) 

The  terrain  factor  rj  is  determined  from  the  following  table  according  to  the  type  of  terrain 
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TABLE  1.  Terrain  Types  and  Terrain  Factors. 


Black 

top 

Dirt 

Light 

brush 

Hard, 

packed 

snow 

Heavy 

brush 

Swampy 

bay 

Loose 

sand 

Soft 

snow 

Terrain  number 

1 

2 

3 

4 

5 

6 

7 

8 

Terrain  factor  rj 

1.0 

1.1 

1.2 

1.3 

1.5 

1.8 

2.1 

1.3  +  0.08 

The  final  equilibrium  heart  rate  HR,  (beats/min)  in  a  given  cycle  is  given  by  (2). 


(  65 

HR,  =<  65  +  0.35  (lhr  -  25) 

M35  +  45  [1-  Exp[-0.01(Ifcr  -  225)]] 

where 

I/,r  =  100(trec.e,/6- 36.76) +  0.4  Wex 


hr  <  25 

25  <  lfcr  <  225 

hr  >  225 


(4) 

(5) 


SWEAT  LOSS  PREDICTION 

The  general  equation  for  predicting  sweat  loss  response  (Am(J  as  a  function  of  exercise, 
environment,  and  clothing  interaction  as  proposed  by  Shapiro  et  al.  (11)  is 

Am.Jgm-ah-1)  =  27.9  Ere,  E^0'455  (6) 

where 

Artiju,  =  change  in  body  weight  from  sweat  loss. 


This  formula  can  be  employed  over  a  wide  range  of  Ere?(50  -  360,  W  nr2)  and  Emax(20  -  525,  W 
nr2),  and  it  is  most  applicable  for  predicting  water  requirements. 


WORK-REST  CYCLE  PREDICTIONS 

The  physical  work-rest  cycle  and  the  time  patterns  of  rectal  temperature  and  heart  rate  are  given 
for  three  different  conditions: 
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(a)  the  time  pattern  for  vesting  subjects  under  various  heat  stress  conditions  referred  to 
as  “resting”  Tre(  (resting  rectal  temperature  at  any  time  t)  and  “vesting”  heart  rate 

HR((r); 

(b)  the  elevation  pattern  for  rectal  temperature  during  physical  work  at  the  given 

climatic  conditions  referred  to  as  “working”  Tret  (rectal  temperature  at  any 
time  t  after  beginning  physical  work)  and  “working”  heart  rate  and 

(c)  the  recovery  rectal  temperature  after  cessation  of  physical  work  referred  to  as 
“recovery”  Tret  (rectal  temperature  at  any  time  t.  after  completion  of  physical 
work)  and  “recovery”  heart  rate  HR^recj. 

The  following  governing  equations  have  been  presented  and  discussed  in  detail  elsewhere  (1). 


temperature 


»..fiS8t  (Tret  )• 

Tret  =  Treo  +  0.1  A  Tre  exp[0.4‘-°’5] 


where 


Tre<  =  rectal  temperature  at  any  time  t 
T reo  =  initial  rectal  temperature  (*C) 

ATre  =  difference  between  the  new  final  equilibrium  rectal  temperature  Tre,  in  the  new 
environment  and  Treo 
t  =  time  (h) 

The  exponential  power  factor,  t  -  0.5,  allows  a  30-min  period  for  the  initial  lag  in  resting  rectal 
temperature  change  when  the  elevation  reaches  about  0.1  of  the  total  change. 

Heart  Rate  at  Rest.  (HR^»). 

HRt(r)  =  65  +  (HR,  -  65)(1  -  e-3 <)  (8) 

where 

HR,  =  equilibrium  heart  rate 
t  =  time  in  hours  (h) 

The  heart  rate  at  rest  in  comfortable  conditions  is  assumed  to  be  65  beats/minute. 


Working  Cycle 


Rectal  Temperature  at  Work. 


Treo  +  0.1  ATre  Exp[0.4‘-°-5  ],  t  <  0 


Treo  +  ATre  [1  -  Expf-?  .  Exp(-0.ir  A Tre  t  )]] 


where 


HR|(w)  =  65  +  (HRy  -  65)[1  -  0.8  Exp  [-(6  -  0.03  (HR,  -  65)t] 


(11) 


where 


^recwork  +  0>5  (Ti 


-  T 

recov  recwor ■ 


fc)'  t  <  ldrec 


Tfer)[l  -Exp{-  a(t  -  tdrec)}],  t  >  tdrec, 


Trecot)  ■-=  Treo  +  (Tre/  -  Treo)[l  -  Exp{-2.0  Exp(-  0.17(Tre/  -  Treo)(t  -  g))}] 

Tretu  r.  rectal  temperature  at  the  beginning  of  decrease  VC);  Treu)  is  not  necessarily  equal 
to  Tre  at  the  end  of  work 

e 

T-Ver  =  equilbiriurn  resting  rectal  temperature  (0) 
a  —  time  constant  of  recovery 
=  1.5  [1  -  ExP{-  1.5  CPcyy}] 

^drec  =  °-25  Exp{-0.5  CPe//} 

CP'j,  =  0.27(^)[P,fc  -  ^PJ  +  (^-)(36  -  T0)  -  1.57 


HR<(rec)  =  HRW  -  (HR^  -  HRr)[l  -  Exp(-  kbt)] 


where 


HR^  =  heart  rate  at  the  end  of  the  work  period 

HRr  =  equilibrium  heart  rate  for  resting  at  the  given  climatic  and  clothing  conditions 


Model:  heat  Network:  6  heat 


Figure  1.  A  Micro  SAINT  Network  diagram  printout. 
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k  =  2  -  O.OICHR^  -  HRr), 
b  =  2.0  +  12[1  -  Expp(~  0.3CPey y)]. 

CLIMATIC  AND  CLOTHING  CONDITIONS 

The  climatic/clothing  conditions  parameters,  im  and  clo,  are  to  be  computed  using 
clo  =  cj  VeC/y 


>m  v  A 
clo  3  ef / 


where 


V#//  =  Vair  +  0.004(M  -  105)  m/s 

and  c i ,  c2,  c3  and  c4  are  taken  from  the  following  tables  for  different  type  of  clothing. 


TABLE  2.  Constant  Parameters  for  Different  Type  of  Clothing. 


Clothing  type 

ci 

c2 

C3 

C4 

closcond 

Shorts 

0.57 

0.3 

1.2 

0.3 

1 

Shorts  and  short-sleeved  shirts 

0.74 

0.28 

0.94 

0.28 

2 

STD  (standard  fatigues) 

0.99 

0.25 

0.75 

0.25 

3 

STD  +  OG  (STD  -f  protective 

1.5 

0.2 

0.51 

0.2 

a 

overgarments  over  the  fatigues) 


MICRO  SAINT  PROGRAMS  FOR  DEEP-BODY  TEMPERATURE,  HEART  RATE, 

AND  SWEAT  LOSS 

The  Micro  SAINT  main  computer  program  that  models  predictive  equations  for  deep-body 
temperature,  heart  rate,  and  sweat  loss  of  clothed  soldiers  performing  physical  work  at  various 
environmental  extremes  is  listed  in  the  appendix.  The  user  is  required  to  supply  the  input  information 
in  the  scenario  section  in  the  program.  The  output  consists  of  the  predictive  variables,  and  it  is  slot  -d 
in  the  snapshot  output  files.  We  list  here  the  input  parameters  for  the  model,  the  output,  the  output 
variables  from  the  model,  arid  an  illustration  of  using  the  developed  software. 

INPUT  PARAMETERS  IN  TIIE  MODEL 

The  input  parameters  for  the  model  are  to  be  entered  in  the  scenario  section  of  the  program  as 
follows: 

g  =  grade  of  terrian  (%) 
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w  =  nude  body  weight  (kg) 

1  =  clothing  and  equipment  weight  (kg) 
vra  =  walking  velocity  (m/s-1) 
ad  =  body  surface  area  (m2) 

ta  =  ambient  temperature  (*C)  in  a  homogenous  environment 

pa  =  saturated  water  vapor  pressure  of  air  at  (mmHg) 

tdb  =  dry  bulb  temperature  (“C) 

tsk  =  average  skin  temperature  (*C) 

psk  =  water  vapor  pressure  at  the  skin  (mmHg) 

vair  =  air  speed  (m/s) 

delt  =  duration  of  time  interval  in  houto  lor  execution  of  problem  variables 

treo  =  initial  rectal  temperature  (*C) 

phia  as  relative  humidity  (%) 

adeff  =  effective  surface  area  for  evaporation  (m2) 

trest  =  time  duration  in  hours  of  the  rest  period 

twork  =  time  duration  in  hours  of  the  work  period 

terrain  =  terrain  number  (Table  1)  determining  rj  as  determined  from  the  terrain  table 
according  to  the  type  of  terrain 
trecovry  =  time  duration  in  hours  of  the  recovery  period 

closcond  =  a  parameter  (Table  2)  det.erming  the  corresponding  values  for  Cj,  c2,  c3,  and  c4 
according  to  the  different  types  of  clothing 

OUTPUT  VARIABLES  FROM  THE  MODEL 

The  following  output  parameters  are  executed  and  stored  in  the  Micro  SAINT  snapshot  files  at 
the  end  of  tasks  2,  3,  4,  5  and  6: 

total  t  =  a  time  variable  (h)  with  duration  from  beginning  of  rest  cycle  to  end  of  recovery 
cycle. 

t  ~  a  time  variable  with  (h)  duration  from  beginning  of  a  corresponding  cycle  to  its  end 
treccqib  =  final  equilibrium  temperature  =  Trey 
temprest  =  equilbrium  temperature  for  rest  cycle 
hrrest  =  equilibrium  heart  rate  for  rest  cycle 

hea'tfr  =  final  equilibrium  heart  rate  for  heat  acclimated  individuals  =  HRy 

swtlosrp  =  sweat-loss  response  (gm~2h-1) 

hrtrest  —  heart  rate  during  rest  cycle  at  time  t 

trecrest  =  rectal  temperature  during  rest  cycle  at  time  t 
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rectemp  =  rectal  temperature  for  rest-work  recovery  cycles  at  total  t 
hrtrate  =  heart  rate  for  rest-work-recovery  cycles  at  total  t 
trecwork  =  rectal  temperature  during  work  cycle  at  time  t 
hrtrwork  =  heart  rate  during  work  cycle  at  time  t 
trecov  =  rectal  temperature  during  recovery  cycle  at  time  t 
hrtrecov  =  heart  rate  during  recovery  at  time  t 

The  following  is  a  listing  of  the  snapshots  for  output  variables  executed  at  the  end  of  tasks  2,  3, 
4,  5,  and  6;  other  parameters  can  be  listed  in  the  snapshots. 

SNAPSHOTS  OF  EXECUTION 


(1)  Trigger: 

End  trigger 

(2)  Task/Network: 

2 

equilib 

(C)  Snapshot  file: 

heateqlb 

Variables  to  store: 

(7)  treceqlb 

(8)  swtlosrp 

(9;  heartfr 

(10)  hrrest 

(11) 

(12) 

(1)  Trigger: 

End  trigger 

(2j  Task/Network: 

3 

rest  temp-heart 

(6)  Snapshot  file: 

restth 

Variables  to  store: 

(7)  t 

(8)  hrtrest 

(9)  trecrest 

(10) 

(11) 

(12) 

(1)  Trigger: 

End  trigger 

(2)  Task/Network: 

4 

work  temp-hear 

(6)  Snapshot  file: 

workth 

Variables  to  store: 

(7)  t 

(8)  hrtrwork 

(9)  trecwork 

(10) 

(11) 

(12) 

(1)  Trigger: 

End  trigger 
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(2)  Task/Network: 

5 

recov  temp-heart 

(6)  Snapshot  file: 

recovth 

Variables  to  store: 

(7)  t 

(8)  hrtrecov 

(9)  trecov 

(10) 

(11) 

(12) 

(1)  Trigger: 

End  trigger 

• 

(2)  Task/ Network: 

6 

data  collect 

(6)  Snapshot  file: 

sumheat 

Variables  to  store: 

(7)  totalt 

(8)  hrtrate 

(9)  rectemp 

(10) 

(11) 

(12) 

ILLUSTRATION  OF  USING  MICRO  SAINT  PROGRAMS  AND  ANALYSIS  OF  OUTPUT  1 

The  following  are  input  parameters  for  an  example  run: 

SIMULATION  SCENARIO 

(1)  Event  Time: 

0.00 

(2)  Expression: 

terrain  =  1; 

closcond  =  1; 

g  =  6; 

w  =  80; 

1  =  8; 

vm  =  1.3; 

ad  =  1.8; 

treo  =  37.2; 

, 

pa  =  20; 

phia  =  .62; 

tdb  =  35; 

tsk  —  36; 
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vair  =  0.515; 
psk  =  44; 
treat  =  1; 
twork  =  1; 
trecovry  =  3; 
ta  =  35; 
adeff  =  .8; 
ddt  =  .1; 

The  terrain  and  the  closcond  numbers  are  to  be  chosen  from  Tables  1  and  2,  respectively,  to 
present  the  environmental  conditions  stated  in  the  tables.  The  main  program  assigns  the  terrain  factor 
T)  and  the  climatic  and  closing  conditions  parameters  automatically. 

The  corresponding  output  results  from  the  snapshot  taken  at  the  end  of  task  2  are: 

“clock”  “treceqlb”  “swtlosrp”  “Heartfr”  “hrrest”  “Trigger:”  “Jobs:” 

0.000000  39.457581  789.714111  158.H3024  83.903366  “End”  “2” 

The  corresponding  output  results  from  the  snapshots  taken  at  the  end  of  tasks  3,  4,  5,  and  6  are 
plotted,  using  Micro  SAINT  utilities,  in  Figs.  2,  3,  4,  and  5,  respectively. 

Figure  2  represents  the  rectal  temperature  and  heart  rate  in  a  rest  cycle  for  1  h  duration.  Figures 
3  and  4  represent  the  temperature  and  heart  rate  during  the  work  and  recovery  cycles,  respectively. 
Figure  5  represents  the  temperature  and  heart  rate  for  the  three  cycles;  it  indicates: 

1)  Rectal  temperature  exhibits  a  slow  decrease  during  the  rest  cycle  while  heart  rate  exhibits 
moderate  increase. 

2)  Rectal  temperature  and  heart  rate  show  large  increases  during  the  working  cycle. 

3)  Heart  rate  reaches  the  resting  equilibrium  condition  at  a  faster  rate  during  the  recovery  cycle 
while  the  rectal  temperature  reaches  the  resting  equilibrium  condition  at  a  slower  rate. 

For  given  environmental  conditions,  one  can  run  the  program  for  a  prescribed  duration  for  the 
rest  and  work  cycles  in  order  to  determine  the  maximum  time  before  the  temperature  and  heart  rate 
reach  their  maximum  allowable  conditions. 


11 


4.  CONCLUSION 


We  developed  a  Micro  SAINT  computer  program  for  a  series  of  predictive  equations  developed  by 
the  U.S.  Army  Research  Institute  of  Environmental  Medicine,  (ARIEM),  for  deep-body  temperature, 
heart  rate,  and  sweat  loss  responses  of  clothed  soldiers  performing  physical  work  at  various 
environmental  extremes.  The  user  can  employ  this  program  by  inserting  input  parameters  in  the 
scenario  section  of  the  Micro  SAINT  program. 

This  computational  module  could  be  used  to  help  reduce  casualties  associated  with  environmental 
heat  extremes  and  to  predict  appropriate  work-rest  cycles  and  water  requirements  equivalent  to  sweat 
loss. 
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APPENDIX 


MICRO  SAINT  NETWORK  AND  PROGRAM  PRINTOUTS 


TASK  NETWORK 


Network  Number:  0 

(1)  Name:  heat  (2)  Type:  Network 

(3)  Upper  Network: 

(4)  Release  Condition:  1; 


(5)  First  Sub-job:  start  initial 

(6)  Sub- jobs  (each  can  be  task  or 

networ] 

Number: 

Name: 

Type: 

start 

initial 

Task 

2 

equilib 

Task 

3 

rest  temp-heart 

Task 

4 

work  temp-heart 

Task 

5 

recov  temp-heart 

Task 

6 

data  collect. 

Task 

7 

end 

Task 

Task  Number:  start 

(1)  Name:  initial  (2)  Type:  Task 

(3)  Upper  Network:  0  heat 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  Deviation:  0; 

(8)  Task/s  Beginning  Effect: 

(9)  Task's  Ending  Effect:  exp=2 . 71 ;t=0 ;totalt=0 ;t3=0 ; tcount=0 ; 

(10)  Decision  Type:  Single  choice 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 


equili  (12)  1; 

(14) 

(16) 

(18) 

(20) 

(22) 

(24) 


Task  Number:  2 

(1)  Name:  equilib  (2)  Type:  Task 

(3)  Upper  Network:  0  heat 

(4)  Release  Condition:  1? 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  Deviation:  0; 

(8)  Task's  Beginning  Effect: 

(9)  Task's  Ending  Effect:  if  terrain  ==  l  then  eta=1.0; 

if  terrain  ==  2  then  eta=l.l; 

if  terrain  --  3  then  eta=1.2; 

if  terrain  ==  4  then  eta=1.3; 

if  terrain  ==  5  then  eta=1.5; 

if  terrain  ==  6  then  eta=1.8; 

if  terrain  ==  7  then  eta=2.1; 

if  terrain  ==  8  then  eta=1.38; 

m=l .  5*w+20*  (w+1) *  (1/w)  A2+eta*  (w+1 )  *  ( l .  5*vm*vm+0 . 35*vm*g)  ; 

mrest=l . 5*w+20* (w+1) * (1/w) A2 ;wexrast=0 ; 

wex=0. 098*g* (w+1) *vm; 

vef f=vair+ . 004* (m-105)  ; 

vef frest=vair+ . 004* (mrest-105)  ; 

if  closcond  ==  1  then  clo=0.57*veff A ( (-1) *0.3)  ; 

if  closcond— 1  then  clorest=0. 57*vef frestA ( (-1) *. 3)  ; 
if  closcond  ==  1  then  imperclo=l . 2*vef f A . 3 ; 

if  closcond==l  then  imclorst=l . 2*vef frest A . 3 ; 
if  closcond  ==  2  then  clo= . 74*vef f A ( (-1) *0 . 28 ) ; 

if  closcond==2  then  clorest=. 74*vef frestA ( (-1) *.28)  ; 
if  closcond  ==  2  then  imperclo=.94*veffA ,28 ; 

if  closcond==2  then  imclorst-. 94*vef frestA . 28 ; 
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if  closeond  =-  3  then  clo=0.99*veffA ( (-1) *.25) ; 

if  closcond==3  then  clorest=. 99*vef f A ( (-1) * . 25) ? 
if  closeond  ==  3  then  imperclo=0 . 75*vef f A . 25 ; 

if  closcond==3  then  iraclorst=. 75*vef frestA . 25 ; 
if  closeond  ==  4  then  clo=l . 5*vef f A ( (-1) * . 2 )  ; 

if  closcond==4  then  clorest=l . 5*vef f restA . 2 ; 
if  closeond  ==  4  then  imperclo= . 51*vef f A . 2 ; 

if  closcond==4  then  imclorst=. 51*vef frestA . 2 ; 
itrest=clorest ;  it=clo; 

im=imperclo*clo ; imrest=imclorst*clorest ; 

hrpc=6. 45*ad* (tdb-tsk)/it ;hrpcrest=6. 45*ad*  (tdb-tsk)  /itrest ; 

dryhtexg=0.0011*hrpc; 

ereq= (m-wex) +hrpc ; 

ereqrest=(mrest-wexrest) ; 

emax=14 . 21*  ( iin/it)  *adef  f  *  (psk-phia*pa)  ; 

emaxrest=14 . 21*  ( imrest/itrest)  *adef  f  *  (psk-phia*pa)  ; 

evaphtex=.8*expA (0.0047* (ereq-emax) ) ; 

evaj.rest=.8*expA  ( .  0047*  (ereqrest-emaxrest) )  ; 

treceqlb=36 . 75+0 . 004* (m-wex) +0 . 0011*hrpc+evaphtex ; 

temprest=36. 75+0. 004* (mrest-wexrest) + . 0011*hrpcrest+evaprest : 

swtlosrp=27.9*ereq*emaxA ( (-1) *0.455) ; 

ihr=100* (treceqlb-36. 75) +. 4*wex; 

if  ihr  <=225  &  ihr  >=  25  then  heartf r=05+0 . 35* ( ihr-25) ; 

if  ihr  >  225  then  heartfr=135+45* (l-expA ( (-. 01) * (ihr-225) ) ) ; 

if  ihr  <  25  then  heartfr=65? 

ihrr est= 1 0  0  * ( tempres t- 3  6 . 7  5 )  + . 4  * wexre  s t ; 

if  ihrrest  <=  225  &  ihrrest>=  25  then  hrrest=65+0 . 35* 

(ihrrest-25) ; 

if  ihrrest  >  225  then  hrrest=135+45* ( l-expA ( (- . 01) * 

(ihrrest-225) ) ) ; 

if  ihrrest  <  25  then  hrrest=65; 

(10)  Decision  Type:  Single  choice 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 


(11)  3 

(13) 

(15) 

(17) 

(19) 

(21) 

(23) 


rest  t  (12)  1; 

(14) 

(16) 

(18) 

(20) 

(22) 

(24) 


Task  Number:  3 

(1)  Name:  rest  temp-heart  (2)  Type:  Task 

(3)  Upper  Network:  0  heat 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  Deviation:  0; 

(8)  Task's  Beginning  Effect: 

(9)  Task's  Ending  Effect: 

trecrest=treo+ (temprest-treo) *0. l*expA (0.4A (t-0.5) ) ; 
hrtrest=65+ (hrrest-65) * (l-expA ( (-1) *3*t) ) ; 

if  t  >=  trest  then  treow=trecrest ;  if  t  >=  trest  then  tcount=0 
rectemp=trecrest ; hrtrat e=hrtrest ; 

(10)  Decision  Type:  Single  choice 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 


(11)  6 
(13) 

(15) 

(17) 

(19) 

(21) 

(23) 


data  c  (12)  1; 

(14) 

(16) 

(18) 

(20) 

(22) 

(24) 
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(2)  Type:  Task 


Task  Number:  4 


(1)  Name:  work  temp-heart 

(3)  Upper  Network:  0  heat 

(4)  Release  Condition:  l; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  Deviation:  0; 

(8)  Task's  Beginning  Effect: 

(9)  Task's  E  (ding  Effect: 
trefs=trecec4xo;hrf=heartfr;t2=(t- (58/m) )  ; 

if  t2  >  0  then  trecwork=treow+ (tref-treow) * (l-expA ( (-1) *2*expA ( (-1) 

. 17* (tref-treow) ) *(t-(58/m) ))); 

hrtrwork=65+(hrf-65)*(l-.8*expA( (-1)  '  ( 6-0 . 03* (hrf-65) 


)  *t) )  ; 

if  t2  0  then  trecwork=  treo+  (temprest-treo)  * .  l*expA  ( .  4A  (t3- .  5) ) 

if  t  >*=  twork  then  tcount=0  ;if  tcount==0  then  tl=twork+delt ? 


if  t  >®  twork  then  trew®trecwork; 


rectemp^trecwor k ; hrtrate*hrt rwork ; 

(10)  Decision  Type:  Single  choice 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 


data  c  (12)  l; 
(14) 

(16) 

(18) 

(20) 

(22) 

(24) 


Task  Number:  5 

(1)  Name:  recov  temp-heart  (2)  Type:  Task 

(3)  Upper  Network:  0  heat 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  Deviation:  0; 

(8)  Task's  Beginning  Effect: 

(9)  Task's  Ending  Effect:  hrr=hrrest; 
k=2- . 01* (hrtrwork-hrr) ; 

cpef f=. 27*imclorst*  (psk-phia*pa) +  ( .  174/clorest.)  *  ( 3  6-ta)  -1.57; 

b=2+12* ( l-expA ( (-1) * . 3*cpef f ) ) ; 

tdrec®. 25*expA ( (-1) *. 5*cpeff) ; 

alpha® 1 . 5* ( l-expA ( (-1) *1 . 5*cpef f ) )  ; 

trer=temprest ; 

if  t  <  tdrec  then  trecov=treo+ ( tref-treo^ * ( l-expA ( (-1) *2*exp- 
(( -1)*. 17* ( tref-treow) )*(tl-( 58/m)) )) ; 
tl=tl+delt ; 

if  t  <  tdrec  then  trecov  =trecwork+ . 5* (trecov-trecwork) ; 

if  t  <  tdrec  |  trecov  >=trew  then  trew=trecov; 

if  t  >=tdrec  then  trecov=trew- (trew-trer) * ( l-expA ( (-1) *alpha* 

(t- tdrec) ) )  ; 

hrtrecov=hrtrwork- (hrtrwork-hrr) * (l-expA ( (-1) *k*b*t) ) ; 
rectemp=trecov ; hrtrate=hrtrecov ; 

(10)  Decision  Type:  Single  choice 

Following  Task/Network:  Probability  Of  Taking 

Number:  Namr'i  This  Path: 


(11)  6 
(13) 

(15) 

(17) 

(19) 

(21) 

(23) 


data  c  (12)  1? 

(14) 

(16) 

(18) 

(20) 

(22) 

(24) 


Task  Number:  6 

(1)  Name:  data  collect. 
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(2)  Type:  Task 


(3)  Upper  Network:  0  heat 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0? 

(7)  Standard  Deviation:  0; 

(8)  Task's  Beginning  Effect: 

(9)  Task's  Ending  Effect:  tcount=tcount+delt ;t=tcount 
if  t3  <=  trest  then  flag=l; 

if  t3  >  trest  &  t3  <=  trest+twork  then  f lag=2 ; 
if  t3  >  trest+twork  &  t3  <=  trest-vtwork+trecovry  then 
if  t3  >  trest+twork+trecovry  then  flag=4; 
t3=t3+delt ; totalt«t3 -delt ; 

(10)  Decision  Type:  Tactical 

Following  Task/Network:  Tactical  Expression: 

Number:  Name: 


rest  t  (12)  flag==l; 
work  t  (14)  flag==2; 
recov  (16)  flag==3; 
end  (18)  flag==4; 

(20) 

(22) 

(24) 


f lag=3  ; 


Task  Number:  7 

(1)  Name:  end  (2)  Type: 

(3)  Upper  Network:  0  heat 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  Deviation:  0; 

(8)  Task's  Beginning  Effect: 

(9)  Task's  Ending  Effect: 

(10)  Decision  Type:  Last  task 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 


(23) 


(12) 

(14) 

(16) 

(18) 

(20) 

(22) 

(24) 


Task 
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